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We numerically examine the self-dual solutions of self-intersecting strings immersed in four dimension. 
By fixing a frame of reference we are able to show that self-intersecting open strings admit world sheet 
topologies that are able to support monopole/anti-monopole structures. Furthermore we identify the 
topologies that can lead to q — q production and annihilation and numerically estimate the production 
rate for different instantons. We find an upper bound on the production rate and use that to estimate 
the string tension /i. By examining the pullback of 75 onto the string world sheet we show that regions of 
chiral symmetry breaking can exist where the string world sheet self-intersects with itself lending credence 
to the relationship of self-intersecting strings and instantons. We are also able to show that for torus 
, knots, the intersection number v = 4(p — q) which validates an earlier conjecture. A program written in 

MAPLE is supplied to exhibit the animations of the torus knots in four dimensions. These animations are 
located at the website http://www-hep.physics.uiowa.edu/~bacus/research.html under the "Animation 
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I. INTRODUCTION 
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Quantum Chromodynamics (QCD) is widely regarded as the correct theory of strong interactions. However, in- 
herently non-perturbative effects such as confinement remain hidden in the theory due to the nature of the strong 
coupling of QCD. In 1973, t'HooftQ proposed the parameter N c of the color gauge group be treated as a free pa- 
rameter, and considered the limit N c — > 00 applied to the expansion of a gauge theory. The resulting expansion gives 
Feynman diagrams which possess the same topology as the quantized dual string model with quarks at the string 
ends. Also, WilsonQ] showed that lattice gauge theory in the strong-coupling approximation is a confining theory due 
Oh, to the formation of color charged strings formed by Faraday flux lines. However this strong coupling approximation 
cannot be made in the continuous theory. Later Makeenko and MigdalQ as well as Gervais and Nevue|| were able 
to establish to first order in the lattice spacing an equivalence between a Wilson loop 
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and solutions to the single string Schroedinger equation, provided that the gauge field satisfied the Yang-Mills equations 
of motion. These indirect probes suggests that the effective action for QCD in the strong coupling regime is a string- 
like theory. Although the precise string theory that governs strongly coupled QCD is unknown we expect this string 
theory to possess certain qualitative features that are characteristic of QCD. Polyakovj7| and Balachandran et al.(|] 
examined the effect of adding a term proportional to the extrinsic curvature to the Nambu-Goto action so that vacua 
effects could be incorporated into the string ansatz. The extrinsic curvature term admits self- intersecting immersions 
as solutions and Nair and Mazur|| relate the self-intersections number to the instanton number. 

In this note we will examine these self-intersecting string immersions while in a fixed frame in four-dimensions to 
see if these configurations possess any realistic topological properties that we expect from an effective action of QCD. 



We will focus our attention on the torus knot solutions of Robertson) 13 in order to be explicit. We will show that 
there exists topologies that supports monopole/anti-monopole fields whose flux is constrained to a finite string (tube) 
and in some cases to strings that are infinite in length. Furthermore we will show that there exists topologies that 
can support pair production and annihilation. By this we mean that the immersed world sheets have remnants that 
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appear as finite segments of strings that emerge and dissolve in a fixed frame. These segments are related to an integer 
q that appears in the self-dual solutions that we consider. We are able to show that q = 111 puts an upper limit 
on the production of fermions. We then examine the pullback of 75 onto the world sheet and show that there can 
be chiral symmetry breaking "bubbles" around the points where the immersion self-intersects with itself. This helps 
strengthen the relationship between intersecting strings on the one hand and QCD instantons on the other. Next we 
verify the conjecture of Robertson fl3f| by showing that for torus knot solutions that the intersection number is given 
by v = 4(p— q). We end by providing a program in MAPLE that annimates the torus knot and torus link solutions. 



II. TOPOLOGICAL SUPPORT FOR QCD PROCESSES FROM SELF-DUAL STRINGS 

One of the questions we ask is what topologies are provided by the string to support the gauge fields and quarks 
in four dimensions. We are interested in the classical features of the self-intersecting instantons since these feature 
should appear as first order quantum corrections to the effective action of QCD in the strong coupling regime. In order 
to appreciate the behavior and support of the strings we need to sit in a frame and examine the topologies that the 
immersed string unfolds. Since the self-dual solutions have non-trivial winding number any features we find related 
to the self-intersection will be protected from quantum fluctuations since winding number is a topological propoerty 
of the string configuration. To proceed we will look at the Nambu-Goto action and include modifications introduced 
by Polyakov. In order to minimize the contributions of crumpled string world sheets to the effective action of QCD, 
Polyakov proposed the action following; 

S = fiJ ^gdadr + ^- J ^gK^K^Jadr (1) 

where K^ a is the extrinsic curvature tensor. However the presence of the extrinsic curvature term not only acts as 
a regulator jHM for bosonic strings but also provides an avenue for the introductions of the topological instanton 
solutions. Furthermore it has been shown that the presence of the extrinsic curvature is necessary for the effective 
action of free fermions projected onto the world sheet |jll|. 

The last summand to the above action can be expressed in an alternate form as; 

Si = — I ^g ab d a t^d b t^dadT (2) 

Si = 7^- I V99 ab ((d a t^ T dj^)(d b t^ T d b t^) ± {d a t^d b t^))dadr (3) 
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where t"" = daX^dbX" (4) 

V9 

As with QCD, the action is minimized by searching for self-dual solutions, viz. d a t IJ ' v = d a t*^ v . WheatcrJT^] showed 
that any surface embedded in four-dimensional space-time which is a complex curve is a solution to the equations of 
motion for the extrinsic curvature term. Robertson|D| showed one explicit example of a self-dual string instanton 
is a (p,q) torus knot. These are examples of superminimal immersions into i? 4 |l4j. Pawelczyk has extended the 
classification of instantons and includes the Euler character of the surface as well ]15| . A torus is parametrized by 
two variables, u and v, representing local coordinates on the surface. The torus knot is the curve on the torus defined 
by the constraint u p + v q = 0. Figure 1 shows an example, the (3,2) torus knot. 
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FIG. 1. The (3,2) torus knot. 
For a (p,q) torus knot instanton solution the string vector is 

X" = mz p ),$t(z p ),$t{-z q ),S(-z g )], 



(5) 



where z — r + ia and 9 and -R are the real and imaginary parts respectively. Note that Robertson has shown somewhat 
more generally that a string vector of the form 

X» = p(F(z))MF(z))MG(z))MG(z))} 



with F{z) and G(z) as any functions analytic in z are also solutions to the equations of motion. These solutions 
automatically satisfy the Euclidean string equations 



d T d T X^ + d a d a X^ = 0, 



as well as the Euclidean constraint equations 



and 



daX^drX^ 1 = 0, 



(d^Xn 2 = (d T X»Y 



(6) 
(7) 
(8) 



To study these solutions, we create computer animated plots of the solutions over a finite time interval. We start 
with a particular (p,q) torus knot solution such as the (3,2) knot. This is a topologically non-trivial solution with 
self-intersection number v = 4. The solution takes the form 

X^{a, t) = [3r 2 cr - a 3 , t 3 - 3ra 2 ,a 2 - r 2 , -2ra] = [x, y, z, t] 

X^ is a mapping from (<r, r) space to (x, y, z, t) space. We choose the last coordinate to be our time parameter and 
invert the the last term and replace r in X^(<r, r) with the solution in terms of X 4 = t and a. (One may chose 
other coordinates for the time such as X 2 where the p = 1 cases correspond to the usual choice where r is the time 



parameter.) In the case of the (3,2) torus knot, the inversion of the last term gives r 



: ; substitution into X^ 



gives 



X»{<T,t) 
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For each value of the time X^ = t we generate a plot of the string. These plots are combined to form the animation. 
The ranges chosen were a = [—1, 1], t = [—1, 1]. The open string is symmetric in our case but we are careful not to 
include characteristics that are due to this symmetry of the solution. We have marked the image of a = 1 and a = — 1 
in order to keep track of the string endpoints on the Reimann sheet.. 
There are several characterisitics that we observed: 
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• As opposed to finite string segments, the majority of sigments are semi-inifinte in length. This is due mainly to 
the inversion of the time parameter. We could have chosen X 2 to be the time component. Then for the small 
sector corresponding to p = l,t would be r. However all other sectors would exhibit the presence of semi-infinite 
segments. These segments may be thought of as flux tubes that carry the monopole and anti-monopole gauge 
fields from anoint out to infinity. One can construct gauge invariant objects that can live on these topologies 



• The strings interact by exchanging flux lines. Whenever two semi-infinite segments, say A and B, touch the 
result is that A will pass its segment above the intersection point onto B while sewing the segment that B had 
above the intersection point onto its lower portion. This exchange of flux lines can also take place with segments 
that are finite. 

• For the case when p is even and q = 2 then the image is that of a single semi-infinite segment that is actually 
the image of the string projected back onto itself. In this case if one had assigned a quark at a = 1 and an 
anti-quarks at a = — 1 then the endpoint of this type of image would correspond to a qq bound state with gauge 
field lines moving up and moving down the flux line that goe s out to spatial infinity effectively canceling the 
contribution of the gauge field in the Wilson line (see Eq. ||il] , jl2"f ) . Such configurations are expected to contribute 
to the < qq > condensate. Examples of this are the (2,2), (4,2), (6,2) and (8,2). 

• One of the most interesting observations is the appearance and disappearance of finite segments of lines. These 
segments can appear at one point, interact with some segment by exchanging flux lines and then disappear 
at another point. The emerging segment can have a qq pair attached to the ends providing topology for pair 
productions and pair annihilation. We will discuss this process in some detail shortly. Nice examples are (4,3), 
(5,3) and (1,3) torus knots. 

• Surfaces that do not intersect p = q, appear as infinite or semi-infinte segments. The multiple images of the 
endpoints can be seen embedded in the string segment itself undergoing pair production and annihilation along 
the flux tube itself. This activity can also be seen in torus knots were p is not relatively primed to q. Examples 
are (6,3), (4,2), and (8,4). 

• Finally we note that the image of a — ±1 does not always map into any endpoint of the image. 



FIG. 2. The (4,2) torus link instanton at time index t=-.9. 
When p(> q) and q and are relatively unprimed, a single semi-infinite segment appears. This is actually the image of two 
semi-infinite strings that mapped on top of each other. 



(see Eq.|l||l|). 
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FIG. 3. The (3,2) torus knot instanton at time index t=0.60 

Figure 3 shows the (3,2) torus knot instanton at time t = 0.60. The ends of the string are labeled with points, and 
the string is colored blue to represent a > and red to represent a < in order to visually distinguish these parts. 
The bounding box containing the plot has dimensions x = —3. .3, y = —3. .3, z = —3. .3. In this frame we see the string 
broken in two halves; where part of the string around a = leaves the bounding box and is mapped out to infinity. 



p=4 q=7 Time Index: -0.75 




FIG. 4. The (4,7) torus knot instanton at time index t=0.60 

Here is an example of a knot with negative intersection number. The vertical strands are semi-infinite in length. 
Here the broken strands are an imaging artifact. Just right of the center one can see the emergence of a finite segment. 
This segment will go on to exchange flux lines with the semi-infinite flux lines and then collapse back into the vacuum. 
These finite strands provide support for pair production and annihilation. 
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p=4 



q=3 Time Index : -2 . 



p=4 q=3 Time Index : -1 . 5 



p=4 q=3 Time Index : -1 . 




p=4 q=7 Time Index : -3 . 



p=4 q=7 Time Index:-2.25 



p=4 q=7 Time Index : -1 . 5 
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Figure 5 shows a sequence of frames of the (4,3) torus knot instanton. Initially, we see two string segments; at 
time t w —1.0, however, a pair of points appear. As time continues to evolve these points develop a string segment 
connecting the two. This finite segment interacts with the pair of strings initially seen, then the two original segments 
combine and disappear. This behavior is reminiscent of the early meson string model, in which a meson was seen as a 
qq pair connected by a string of gluon flux. In this interpretation, it appears as though a qq pair is created out of the 
vacuum, interacts strongly with a pair of gluon flux lines, then is annihilated back to the vacuum in the presence of 
a monopole/anti-monopole pair. Here it is the original string segment that we interpret as monopole-anti-monopolc 
pairs. We emphasize that these are open string solutions that provide this topology as opposed to closed stings. The 
new requirement is that the opens strings be self intersecting. 

The structure of the images suggests that the quarks interact with each other through non-local interactions. We 
would like for the strings to correspond to flux tubes of gauge fields and to have quarks on the cndpoints. Lets fix 
the gauge, Aq = and consider the Wilson line integral given by the following path-ordered exponential, 

C a+ dX i 

U(X(a + ,t) > X(a.,t))=Vexp( A,— da). (9) 

J a- d(J 

Here X(a+) and X{a-) are the end points of the string as seen from a fixed frame and a is used to parameterize 
the string for each time, t. From the image of the string immersions it appears that non-local point particle-like 
interaction Lagrangians such 

L = q(X(a + ,t))U(X(a + ,t),X(a.,t))q(X(a^t)). (10) 

can be supported. In order to account for the semi-infinite segments we can imagine taking one of the endpoints of 
the Wilson line out to spatial infinity. At spatial infinity the gauge fields approach zero and the Wilson line goes to 
the identity. The semi-infinite segment then can support a quark attached to the end of the segment with a monopole 
or anti-monopole gauge field confined to the segment. One can build an operator that respect the residual gauge 
invariance with 

q(X(a + ,t))U(X(a+,t),<x>), (11) 

or either 

U(oo,X(a.,t))q(X(a^,t)). (12) 

The lines that are infinite in extent correspond to both Wilson line endpoints being taken to spatial infinity. These 
configurations seem to be the catalyst for pair production and annihilation. There might be suppression to these 
infinite and semi- infinite segments if one uses the same "area law" arguments j^j one uses to suppress largely separated 
quark paths in pair production if the flux tubes are too far from each other. However very close pairs q U and U q 
may be significant and could contribute to < qq >. To sum up the overall picture that these string suggests in that 
the general vacuum structure will consist of infinite and semi-infinite flux tubes as well as emerging and disappearing 
finite tubes and stable finite strings. The gauge fields will be confined to the surfaces of these segments suggesting 
that monopole and anti-monopole are pervasive and where pair production seems to be catalyzed by the presence of 
these monopoles/anti-monopoles. The presence of the flux tubes pouring their gauge flux out at spatial infinity is 
already suggested by the electric Meissner effect Jl6| . The exchange of flux tubes between pair-produced segments and 
the other segments suggest the charge-exchange scattering processes between fermions and charged Dirac magnetic 
monopoles jf7[. 



III. QQ PRODUCTION 

The topology that supports the pair production and annihilation is not due directly to string self-intersection but 
is related to how the time parameter X 4 takes its values from the a — r parameter space. This is related to the fact 
that the solution is self-dual, however. Let us consider the multi- valued relationship of r with a for a given time, X±. 
Since the X4 coordinate only depends on g, we can examine a time sequence of the relationship of a vs. r. To be 
explicit let us consider any solution where q — 3. 
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FIG. 7. a vs. r at different times for q = 3 

Here one sees that at very early times (t << — 1), t is a single valued function for values of a between —1 and 
1. As time moves on, this function becomes multi- valued (here at t = —2) for values where a takes values from 
-1 to 1. Now the physical a — r parameter space is split into three distinct regions corresponding to three distinct 
images in the frame. As time moves on these region move toward the origin until at t = all three regions merge 
for an "instant" and the regions begin to recede. At some later time (here t — 2) the r parameter is again single 
valued with respect to the physically relevant values for a. This picture highlights two issues about the self-dual 
solution that may be realized in the strong coupling regime of QCD. The first is how the multivalued nature of r vs. 
<t can cause pair production and annihilation and the second being the singularity of the map at t = 0, the point of 
self-intersection, leads to regions where chiral symmetry can be broken. These singular points are where the string 
segments are allowed to exhange flux lines and also where the determinant of the induced metric vanishes. The details 
of this last remark will be explored in the next section. For now we will take a closer look at the issues related to 
production and annihilation. 

From figure 0, one sees that from the time of the production (t = —2) until the time of annihilation (t = 2), four 
units of time have transpired. Here we can assume that the constant /i from the action in Eq. [1], sets the space and 
time scale since fi has dimensions of [£~ 2 ]. From the uncertainty principle we can estimate a typical value of 

since four quarks are produced in a time 4;4=. This would imply that /i w (16M g ) 2 where M q is the quark mass. The 
q = 3 set of torus knots correspond to the minimal pair-producing configuration. However higher q instantons could 
easily surpass the bounds from the uncertainty principle by producing more quarks in nearly equal times to that of 
the q — 3 solutions. For this reason we need to check the production capabilities of the q > 3 solutions. 
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In order to determine the time it takes for higher q configurations to produce and annihilate quarks we need only 
to ask at what time, t rea i are all the roots of t(g = — 1, t rea i) are real. For every value of q there are q regions that 
must merge into the (— 1 < a < 1) region and then move out again. The total production/annihilation time would 
then be 2 x t rea i. As an example lets examine the q = 4 case. 



-1.5= Time -0.46= Time -0.2= Time 




FIG. 8. a vs. r at different times 

In this case the pair production begins at t = where r takes values of and -^=. The time for total 

pair production is then time units or about three units of time. In this amount of time again four quarks or 
anti-quarks are produced and annihilated. As q increases it empirically appears that the amount of time to undergo 
full pair production monotonically decreases but reaches a plateau. One can show that the total production times for 
different values of q are: 



5 = 3, 


treai 2 


(13) 


5 = 4, 


+ 8 

Weal - 3^ 


(14) 


q = 5, 


^real = 1.37 


(15) 


q = 6, 


Ueal = 1-285 


(16) 
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q = W, t real = l.U7 (17) 

q = 100, t reai = 1.012 (18) 

9 = HI, t rea i = .608 (19) 

We observe that 9=111 puts an upper limit on the value of q that will produce any new pairs. For all values greater 
that q = 111 some of the roots for r remain complex even when t = 0. By analogy with the q = 3,4 and 5 cases 
we see that there can be at most 112 quarks/anti-quarks produced in 1.2 seconds by any instanton. Since the total 
number of quarks is limited by all torus knot instantons we can again use the uncertainty principle to get an even 
better estimate of fi. We find that /Lt w (112M 9 ) 2 . 

A nimations of several (p, q) solutions wi th AT 4 chosen as the time parameter can be found at http://ww 



hep . physics . uiowa.edu / ~ bacus/research . html . 



IV. CHIRAL SYMMETRY BREAKING PROCESSES 

As we discussed earlier, the self intersection of the world sheet can cause singular points to exist and lead to zeros 
in the determinant of the induced metric g a ^ = d a X p d\ } X u . Self duality for the string instantons imply that t^ u from 
Eq.([|) satisfy the constraint 

= ^ vXp txp + C tlv , (20) 

where C 1 ' 2 = — C 2,1 = C 3,4 = — C 4 ' 3 = 1 and all other entries are zero. 

As we are interested in chiral symmetry, we would like to look at the pullback of 75 onto the string world sheet. In 
terms of functions on the world sheet we may write 

75 = ^e^d^d^N^^^} [ 7A , 7p ]. (21) 

In the above N% are elements of the normal bundle of the string world sheet and e AB is the volume form associated 
with the normal bundle. By using self duality we can eliminate the tangent vectors from the above expression and 
write, 

75 - h EF C^N x E N p F e^ Xplb . (22) 

The normal vectors satisfy the conditions that d a X p N E = 0, and N^N flF = 5^,. These conditions as well as the 
self-duality conditions imply that a solution for the tangent and normal vectors is 

d a x» = (d a x\d a x 2 ,d a x 3 ,d a x 4 ) (23) 

d T X» = (-d a X 2 ,d a X 1 ,d a X*,-d a X 3 ) (24) 
iVf =- t =(a CT X 3 ,a CT X 4 ,-a CT X 1 ,-a CT X 2 ) (25) 
N$ = -^(dvX^-dvX^doX^-drX 1 ), (26) 

where g is the determinant of the induced metric. For these solutions *Jg = g aa . Without loss of generality we can 
take the (3,2) torus knot as an example. Here the torus knot determinant is 

= 9(t 4 + a 4 ) - 12tV + 4<t 2 + 4r 2 . 

Clearly this vanishes at a = 0,t = 0. This corresponds to the point where the string self-intersects. At this point 
the normal bundle is ill-defined and two copies of the tangent vectors are imaged at this point. The pullback of 75 is 
also ill-defined when the determinant vanishes. Elsewhere on the manifold 75 can be used to construct a projection 
operator for chiral symmetry. However these points where the determinant vanish act as tiny "bubbles" where chiral 
symmetry is broken. The transition of the normal bundle going through this singularity has the effect of imposing 
reflections on some but not necessarily all of the components. In the (3,2) case the relationship of the normal vectors 
just before and just after is given by: 

N? bef ore = -^(^ 3 ,^V^A^-d ff X 2 ) (27) 
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while 



N? after = — {d«X\-d a X\-d a X\d a X 2 ) (28) 

flipping the second and fourth components. Similar changes happen in both the tangent vectors and the other 
normal. To show that the orientation has changed consider the vector M = Ni x N2 for the (3,2) knot. Just after the 
singularity, M y — > — M y while all other components remain the same. This corresponds to a change in the orientation. 
This pictures corresponds to the field theoretic case where the core of monopoles can be thought of as bubbles where 
chiral symmetry is broken. The S-wave of the fcrmions interact with this core which leads to processes such as those 
seen in 0,|§. 



V. THE SELF-INTERSECTION NUMBER 

It is conjectured in Jl3[ that the self-intersection number of the (p,q) torus is v = 4(p — q). Here we explicitly 
calculate the self-intersection number v for a general (p,q) torus-knot solution proving the conjecture. Consider an 
open string where (-L < a < L) and (—00 < t < 00). Starting from the general definition of the self-intersection 
number 

-1 /*oo pL 

v = — / dr dadj^dat^ 



2tt 



— OO 



— L 



which we can express in terms of the string solution X(cr, r) as 

1 poo pL cd f c.d 

v = ^r\ dT dad a (—d c x»d d x»)d a (—d c x»d d xn 

Using the general (p,q) torus-knot solution 

X(a, t) = [(a 2 + T 2 f' 2 sin(ptan- 1 (-)), (a 2 + t 2 )p' 2 cosfptan"^-)), 

T T 

-(a 2 + r 2 f/ 2 cos(q tan-^-)), -{a 2 + T 2 f' 2 sin^tan"^-))] 



we find 



1 f°° , f L , 8p 2 q 2 {p-q) 2 (<7 2 +t 2 )p+i- 1 
v = — / dr d<7- 



If we let r 2 = a 2 + r 2 



v = — 



2W-oo J-L (p 2 (<T 2 +T 2 )P + q 2 (a 2 +T 2 )i) 2 



1 /-Va f^koj g p 2 q 2(p^ q )2 r 2(p+ q -l) + l 

' dO I dr- 



k J-tt/2 Jo (p 2 r 2 P + q 2 r 2 i) 



2 



,/a (p~q)q 2 (q 2 r 2 (P~i) +p 2 ) 



coa(8) 




Let 6 = 29 + 7r, 



to) r to ^-2(p- ff) 

^ 7o (£X(p-g)) 2 ( 1 ~ c ° sW )(p-g) + 1 



Note the following relation; 

1-2 r. 

/Ysin Q. cos 0)d0 = -i J f(— 

2i 

Then 



/(sin0,cos0)d0 = -i j> /(- 



12 



v = ^ p - q) I * , 2 (p-g) 

= 4(g — p) — — ^ — — 7r Residues in unit circle 
= 4(g — p)(l — Residues in unit circle) 

Now we must find the singular points of f(z) = — — rr ; they are 

JK ' *((fi"-«) 2 (i(i- s± §- ! -)''-«+i)' 

z = 



z = l- 2(-(^) 2 ( J L) 2 ^ <? ))^ ± 2./(-(^) 2 (L) 2 (P-2))i^ + (i(L)(P"9))^F 
P V V V 

Evaluating the residue [z — Zo)f(z)\ z — Zo gives zero; then we have 

v = i(q-p) 

for (p,q) torus-knot solutions. This confirms the conjecture of Robertson fll3| |. 

VI. COMPUTER PROGRAM 

The software described in this section was used to implement and generate the animated solutions and to perform 
some lengthy calculations. Two symbolic computation systems were used; Maple and Mathematica. 

Listing 1 shows the Maple source code which generates raw (x, y, z) coordinates for a (p,q) torus knot solution and 
exports them to a file. The user selects the values for p and q, the frame numbers to start and end at, the number of 
data points to use for each frame, and the ranges for the parameters t and a. 



Listing 1 



# Animation of the string instanton solutions. This program generates a# 

# .dat numeric data file for each frame of the animation sequence. # 

# Usage (on UNIX): maple 4D_TorusKnot_Data.txt & # 

# # 

# Bob Bacus c. 1997 # 

# User input parameters # 

p : =3 : # F=zp # 
q:=2: # G=-zq # 

startf rame : =0 : # Frame to start from # 

endf rame : =400 : # Frame to end on # 

points: =400: # Number of data points per frame # 

lmin:=-3:lmax:=3: # Range for the time parameter # 

smin: =-1 : smax: =1 : # Range for the length parameter # 

# Main Program # 

gc (300000) : interface (screenwidth=500) : Digits : =4 : 
readlib(unassign) :readlib (write) : 

unassign( ' t ' ) : unassign( ' 1 ' ) :unassign( 'n' ) :unassign( 'u' ) :unassign( ' s ' ) : 

assume (t, real) : assume (s, real) : 

z:=t+I*s: 

F:=expand(zp) : G : =expand(-zq) : 

imf :=Im(F) :ref :=Re(F) :reg:=Re(G) :img:=Im(G) : # X=[Im(F) ,Re(F) ,Re (G) , Im(G) ] 
l:=(lmax-lmin)/ (endf rame-startf rame) *n+lmin: # n represents frame number; 1 is the time 
v:=(smax-smin)/(points-l) : 
for n from startframe by 1 to endframe do 

open (cat (convert (p, string) ,'x', convert (q, string) , ' , convert (n, string) , ' . dat 1 ) ) : 

writeln(evalf (1) ) : 

writeln(points) : 

for i from 1 to q-1 do 

writeln(cat ( 1 Solution #' , convert (i , string) ) ) : 
for u from smin by v to smax do 

r:=[fsolve(subs(s=u,img)=l,t, complex)] :r:=r[i] : 

X:=expand(subs (t=r , imf ) ) : Y: =expand(subs (t=r ,ref ) ) : Z :=expand(subs (t=r ,reg) ) : 
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writeln(evalf (subs(s=u, [X,Y,Z] ))) : 
od; 

od; 

closeO : 

od; 

unassign( 't ' ) : unassign( ' 1 ' ) :unassign( 'n' ) :unassign( 'u' ) : 
closeO : 



Listing 2 

# 

# Animation of the string instanton solutions. This program generates a 

# .gif graphics file for each frame of the animation sequence. 
# 

# Bob Bacus c.1997 

# User input parameters 

p:=l: # F= zp 
q:=2: # G=-zq 

st art frame : =0 : # frame to start with 
endf rame : =400 : # frame to end with 

# 

space: =< ' :gc(300000) : interf ace(plotdevice=gif ) : 
readlib (write) :readlib(unassign) : 
for n from startframe by 1 to endframe do 
d: = [] : 

file:=cat(convert(p, string) , 'x' , convert (q, string) , , 

convert (n , string) , 1 . dat ' ) : 
indx:=parse(readline(f ile)) : 

if indx>=0 then tp : =cat ( 'p= ', convert (p, string) , ' q= ', convert (q, string) , ' Time Index : 
' , substring ( convert (indx, string) , 1 . . 7) , substring (space , 1 . . 7-length (substring (convert ( 
indx, string) , 1 .. 7) )) ) else tp:=cat( 'p=' .convert (p, string) , ' q= ', convert (q, string) , ' 
Time Index : ' , substring (convert (indx, string) , 1 . . 8) , substring (space , 1 . . 8-length( substring ( 
convert (indx, string) , 1 . . 8) ) ) ) f i : 
points : =parse (readline (file) ) : 
for i from 1 to q-1 do 
readline (file) : 
j:=l: 

for s from 1 to points/2 do 
pt :=parse (readline (f ile)) : 

if s=l then if linalg [iszero] (map ( Im.pt)) then 
smin:=pt :b(j):=[pt] elseb(j):=[] : smin:=[2.999,2.999,2.999] fifi: 

if s>l then if linalg [iszero] (map(Im.pt)) then b(j) :=[op(b(j)) ,pt] else j:=j+l :b 
(j): = [] : fi fi: 
od; 
k:=l: 
c(l): = []: 

for s from points/2+1 to points do 
pt : =parse (readline(f ile) ) : 

if s<points then if linalg [iszero] (map(Im.pt)) then c(k) :=[op(c(k)) ,pt] 
else k:=k+l: c(k):=[]: fifi: 

if s=points then if linalg [iszero] (map(Im.pt)) then smax:=pt: c(k) :=[op(c(k)) ,pt] 
else smax:=[2.999,2.999,2.999] fi fi: 
od; 

Ll: = [] : 
L2: = [] : 

for h from 1 to j do 

if nops(b(h))>0 then LI : = [op (LI) ,b(h)] fi: 

od: 

for h from 1 to k do 

if nops(c(h))>0 then L2:=[op(L2) ,c(h)] fi: 

od: 

if nops(Ll)>0 then LI : =op(l ,L1) : else Ll:=[]: fi: 

if nops(L2)>0 then L2:=op(l,L2) : else Ll:=[]: fi: 

d:=[op(d) , PL0T3D (POINTS ( [2.999,2.999,2.999] ) ,CDL0R(RGB, 0,0,0)) , 
PL0T3D (POINTS ( [smin , smax] ) , C0L0R(RGB ,0,0,0)), PL0T3D (CURVES (LI) , COLOR (RGB ,1,0,0)), PL0T3D (CURVES (L2) , 
C0L0R(RG B, 0,0,1))] : 
od: 

interface (plotoutput=cat (convert (p, string) , 'x' , convert (q, string) , 
'-' .convert (n, string ),'.gif')): 
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plots [display] (d, view= [-3 . . 3, -3 . . 3, -3 . .3] ,projection= . 1 .orient at ion= [25, 60] ,axes=boxe 
d,labels=['x' , 'y' , 'z'] ,scaling=constrained,titlefont=[CDURIER,10] ,title=tp) ; 
test : =readline (file) : 

unassign( ' 1 ' ) : unassign( 'b' ) :unassign( 'points ' ) : close () : 
od: 
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